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' Abstract 

I In this paper, we consider the inverse problem of recovering a doubly periodic Lipschitz 

• structure through the measurement of the scattered field above the structure produced by 

I point sources lying above the structure. The medium above the structure is assumed to 

be homogenous and lossless with a positive dielectric coefficient. Below the structure is a 
perfect conductor partially coated with a dielectric. A periodic version of the linear sampling 
method is developed to reconstruct the doubly periodic structure using the near field data. 
In this case, the far field equation defined on the unit ball of M.^ is replaced by the near field 
equation which is a linear integral equation of the first kind defined on a plane above the 
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' periodic surface. 
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1 Introduction 

^ , Consider the problem of scattering of electromagnetic waves by a doubly periodic structure of 

$H ' period A = (2ir, 2ir) defined by 

r = {x3 = /(xi, X2) I f{xi + 2ni7r, X2 + 2n27r) = /(xi, X2) > 0, V n = (ni, 722) G Z^}, 

where the function / is assumed to be Lipschitz continuous so the periodic structure F is a 
Lipschitz surface. The medium above the structure is assumed to be homogenous with a constant 
dielectric coefhcient eo > 0, and below the structure is a perfect conductor with a partially 
coated dielectric boundary. The magnetic permeability is assumed to be a positive constant fiQ 
throughout M^. Given the structure and a time-harmonic electromagnetic wave incident on the 
structure, the direct scattering problem is to compute the electric and magnetic distributions 
away from the structure. In this paper, we are interested in the inverse problem of reconstructing 
the shape of the bi-periodic structure from a knowledge of the incident and scattered fields. The 
purpose of this paper is to develop a periodic version of the Linear Sampling Method for such 
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an inverse problem. We refer to [24] for historical remarks and details of the applications of 
the scattering theory in periodic structures and [8] for a recent overview of the linear sampling 
method. 

Physically, the propagation of time-harmonic electromagnetic waves (with the time variation 
of the form e~*'^*, w > 0) in a homogeneous isotropic medium in is modeled by the time- 
harmonic Maxwell equations: 

cmlE -ikH = 0, cuvlH + ikE = 0. (1.1) 

Here, we assume that the medium above the structure is lossless, that is, A: is a positive wave 
number given by A; = y^ioT^ow in terms of the frequency w, the electric permittivity Eq and the 
magnetic permeability fiQ. Consider the time-harmonic plane wave 

incident on T from the top region := {x G | 2:3 > f{xi,X2)}, where d = {ai, 02, —(3) = 
(cos 9i cos 62, cos 9i sin 62, — sin ^i) is the incident direction specified by 9i and 82 with < < 
vr, < ^2 ^ 27r and the vectors p and q are polarization directions satisfying that p = y^jL/e{q x d) 
and qJ-d. 

In this paper, we assume that the boundary F has a Lipschitz dissection T = r/jUSUF/, 
where Td and Tj are disjoint, relatively open subsets of T having S as their common boundary. 
Suppose below F is a perfect conductor partially coated by a dielectric on Tj. The problem of 
scattering of time-harmonic electromagnetic waves is modeled by the following exterior mixed 
boundary value problem: 

curl cvltIE - k'^E = in Q, (1.2) 
uxE = on Td, (1.3) 
u X cuilE - iX{i^ X E) X I' = on Tj, (1.4) 
E = E' + E' in VL, (1.5) 

where v is the unit normal of F pointing into VL. We assume throughout this paper that A is a 
positive constant. 

Set a = (ai,a2)0) G M^, n = (ni,n2) G Z^. Motivated by the periodicity of the medium 
we look for a-quasi-periodic solutions in the sense that £'(xi, X2, X3)e~*"'^ is 2tt periodic with 
respect to xi and X2, respectively. Since the domain is unbounded in the xs-direction, a radiation 
condition must be imposed. Physically it is required that the scattered fields remain bounded as 
X3 tends to -|-oo, which leads to the so-called outgoing wave condition of the form: 



^ ^„e^(""-^+^"^3)^ 3:3 > max/(xi,X2), (1.6) 

neZ2 




where a„ = (ai + ni, 02 + "2, 0) G M?, E^ = (E^^K Ei?\ E^n^ ) G are constant vectors and 

if I On I < k, 
if |Q!n| > k, 

with i'^ = —1. Furthermore, we assume that /?„ / for all n G Z^. The series expansion in (II. 6p 
will be considered as the Rayleigh series of the scattered field, and the condition is called the 
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Rayleigh expansion radiation condition. The coefficients En in (jl.6p are also called the Rayleigh 
sequence. From the fact that divE'^(x) = in 17 it is clear that 

an-En + (3nE(^^ =0. (1.7) 

The inverse problem considered in this paper is concerned with determining the profile F 
and the impedance coefficient A from a knowledge of the incident wave E^ and the tangential 
component of the total electric field, i' x E, on a, plane F^ = {x G | X3 = 5} above the 
structure. The uniqueness of this inverse problem was proved in |15j for the case when the incident 
waves are electric dipoles. Precisely, it was shown in |15) that, if the tangential components 
on F;, of two scattered electric fields are identical for all incident electric dipoles S*"'(x;y) = 
curla;Curla;{PG(a:, y)} with all y G F^ and three linear independent vectors P G M^, then their 
corresponding scattered periodic structures Tj{j = 1,2) and the impedance coefficients Xj{j = 
1,2) on Tj must coincide, where G{x,y) is the free-space quasi-periodic Green function (see 
Section[2]). In this paper, we are interested in numerically reconstructing the shape of the periodic 
structure F by using the idea of the linear sampling method. The linear sampling method was 
proposed in [l3] for numerically reconstructing the shape and location of the obstacle in the 
inverse acoustic obstacle scattering problems. This method has attracted extensive attention in 
recent years since it does not need to know the physical property of the scattering obstacles in 
advance. The application of the linear sampling method to the inverse electromagnetic scattering 
problems can be found in [9| I10|, [T2]. Recently in [I6l, a periodic version of the linear sampling 
method was proposed and implemented for the two-dimensional TE polarization case of the 
inverse problem considered in this paper, where the Maxwell equations are replaced by the scalar 
Helmholtz equation and the boundary conditions on F^) and Tj are replaced with the Dirichlet 
and impedance conditions, respectively. In |17j . Kirsch proposed a mathematically-justified 
version of the linear sampling method, the so-called factorization method. However, it is still an 
open question to characterize a bounded conducting obstacle for the Maxwell equations by using 
the factorization method (see [18]). We refer to [2][3l[l9] for the application of the factorization 
method to the 2D inverse problems by diffraction gratings with the Dirichlet, impedance and 
transmission conditions and to [1] for a recent convergence result of the linear sampling method 
as well as a connection between the linear sampling and factorization methods. 

The inverse scattering problem by a smooth doubly periodic structure has been studied in 
[U [5] for the case when F/ = 0. With a lossy medium (i.e., Im (fe) > 0) above the conductor, 
Ammari [1] proved a global uniqueness result for the inverse problem with one incident plane 
wave. For the case of lossless medium (i.e., Im (A;) = 0) above the conductor, a local uniqueness 
result was obtained by Bao and Zhou in [5] for the inverse problem with one incident plane wave 
by establishing a lower bound of the first eigenvalue of the curl curl operator with the boundary 
condition (|1.3p in a bounded, smooth convex domain in M^. The stability of the inverse problem 
was also studied in [5]. Recently in [6] it was proved that one incident plane wave is enough to 
uniquely determine a bi-periodic polyhedral structure except for several extremely exceptional 
cases. 

Note that the inverse problem we are concerned in this paper involves in the near field 
measurements since only a finite number of terms in (jl.6p are upwards propagating plane waves 
and the rest are evanescent modes that decay exponentially with distance away from the grating. 
Thus we use near field data rather than far field data to reconstruct the grating structure, which 
implies that the far field equation defined on the unit ball of for the non-periodic case must 
be replaced by a near field equation defined on a plane above the structure. On the other 
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hand, instead of using electromagnetic Herglotz pairs in the case of bounded obstacle scattering 
problems, we consider another kind of incident electric fields (see Section [3] and Remark 14. 2p 
which lead to a denseness range result on the grating structure since scattering occurs in a 
half space and the solution is a-quasi-periodic depending on the incident angle of the incident 
direction. This differs from the original version of the linear sampling method which makes use 
of incident plane waves of all incident directions with three linearly independent polarization 
directions (cf. [9j). 

The remaining part of the paper is organized as follows. Section [2] is devoted to the basic 
quasi-periodic function spaces used in the study of electromagnetic scattering problems by peri- 
odic structures. Section [3] gives several important lemmas which are necessary for establishing 
the main result. The main result on the periodic Linear Sampling Method and the numerical 
strategies on the implementation of the linear sampling method are presented in Section HI 

2 Basic Function Spaces 

In this section we introduce some quasi-periodic Sobolev spaces which are well-suited for our 
problems. Due to the periodicity of the problem, the original problem can be reduced to a 
problem in a single periodic cell of the grating profile. To this end and for the subsequent 
analysis, we reformulate the following notations: 



for any b > max{/(xi, a;2)}. We now introduce the scalar quasi-periodic Sobolev space: 




n 



r 



{x-i = f{xi,X2) I < Xl, X2 < 27r}, 

{x3 = 6 I < xi, < 2tt}, 

{x G I X3 > /(xi, X2), < xi, X2 < 2tt} 

{x G I X3 < 6}, 

{x G : < xi,X2 < 27r} 



= {u{x) = exp[i(a„ • x + /3„X3)] | u E L'^{nb),Vu £ {L^{VLi,)f ,Un € C}. 



nGZ2 



Denote by if 2 (r^) the trace space of H^{Q.i)) on F;, with the norm 



n ^ tl hi 



where /„ = (/, exp(ia„ • x))i2(Y^y Write H 2 (r^) = {H-2{Tti)y, the dual space to H^{Tb). We 
also need some vector spaces. Let 



iJ(curl,J7fc) = {E{x)= ^ Enexp[i{an-x + /3nX3)]\En£C^, 



E e curlE e {l?iS^b)f} 



with the norm 





+ 1 1 curl £^1 
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and let 

i/o (curl, Oft) = {E e H {curl, fib), x = on TJ. 

Define 

X := X{nb,ri) = {EG H{cml,nh), V X ^Ir, G ^^(r,)} 

with the norm 

ll^llx = ll^llH(curl,n,) + ll'^xS||2,(r,)' 

where L\{T) = {E ^ {L?{Y:)f, v-E = OonT}. For x' = (xi,X2,6) G Tfc, s G M define 
Hl^Tb) = {Eix') =Y.^n exp(ia„ • x') | ^„ G C^ 63 • ^ = 0, 

nGZ2 

nGZ2 

i/^div , Tfe) = {E{x') =Y,En exp{ian • x') | K G C^ 63 • ^ = 0, 

ll^lln^div.r,) = E (1 + + \En ■ a„p) < +00} 

nGZ2 

/7,^(curl , Tfe) = {E{x') =J2^^ exp(ia„ • x') | G C^, 63 • = 0, 

nGZ2 

ll^llH^curlA) = E (1 + + 1^- >< < 

neZ2 

and write (Fi,) = i^°(rfe). Recall that 

F-^/'(div,r,) = {63 X E\rjEe H{cml,nt)} 

— 1/2 

and that the trace mapping from //(curl , ilt) to (div , Fh) is continuous and surjective (see 

[7] and the references there). The trace space on the complementary part To of X{il.f),^l) is 

Y(Td) = {/ G {H-^/\rD)f\3E G (curl, r?b) such that X E|r, G L2(r,), z. x ^|r^ = /} 

which is a Banach space with the norm 

ll/lly(r,) = ^^mE\\jj^curi,n,) + \W^E\\l,ir,)\ 

E G Hoicml, fib), y X ^Ir, G L\{Vi), v x ^|r^ = /}. 

An equivalent norm to || • ||yp^ is given by (see [HI [HI 122) ) 

lll/llli= sup ii^^if'^ \ 

where, for E G f/o(curl , 0^) satisfying that v x i?|r/ G L^(r/) and x -EItj, = /) we have 

<f,V>i=[ {cmlE-V-E-cmlV)dx- [ uxE-Vds, e X{flb,Ti). (2.1) 
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In particular, Y{T£)) is a Hilbert space, and (j2.ip can be considered as a duality between Y(T£)) 
and its dual space y(r£))'. From (I2.ip it can be seen that (/? G y(r£))' can be extended as a 
function G -^curi'^(-'^) defined on the whole boundary T such that ^Ir^ € Lf{Tj). 

We conclude this section with introducing the following free space Green function ^>(x, y) for 
the Helmholtz equation (A + k'^)u = in M^: 

^ik\x-y\ 

'^{x,y) = J-. 1 

4:7r\x — y\ 

and the following free space a-quasi-periodic Green function G{x, y) for the Helmholtz equation: 
y) = ^ y2 exp(ia„ ■ {x - y) + il3n\xs - ys\) (2.2) 

neN2 ^ 

with an,f^n defined as in the introduction. 



3 Several Lemmas 



In this section we prove several important lemmas which are necessary for the proof of the main 
theorem. We first define the incident electric field 



S™(x;5) := curl^curl^ / 9iy)G{y,x)ds{y) 



(3.1) 



for g G L^(rb). From the definition of G{x,y) it is seen that E^^(x]g) satisfies the radiation 
condition (II. 6p in the region below Ffe. This means that physically the above incident field 
propagates upward and does not appear to be meaningful as incident waves. Thus the total 
electric field corresponding to E^^{x;g) can not be generated directly. We will discuss how to 
solve the direct scattering problem for such incident waves in the final section. For any g G L^(Ff,) 
we next define a function (Hg) G B := Y{T£)) x L^(F/) by 



u{x) X E'''^{x;g), on Fd, 

^{x) X cur\E^"{x; g) — iXE^"{x; g)T, on Fi, 



where, for any vector field V , Vr ■= {i^ x V) x u denotes its tangential component on a surface. 
Lemma 3.1. The range of H is dense in B. 

Proof. For f x h £ B* := Y{T£)) x L^(F/), we are going to prove that / = 0, /i = under the 
assumption that < Hg,f x h >b,b* = for any g G L^(Ff,). Recalling that the duality between 
Y{Td) and Y{rn)' is defined by dHJ and the duality between L'f{Ti) and L^(F/) is the 
scalar product, we have 

• f{x)ds{x) 



/ iy{x) X 


curl curl 2. / 




Jv 



+ j \ v{x) X 



curUcurUcurU / g{y)G{y,x)ds{y) 



+iX 



curUcurU / g{y)G{y,x)ds{y) 



} • h{x)ds{x) 
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Since / G Y{T£)) , there is an extension / G -f^curl °^ defined on F, satisfying that 
/Itj G -^((r/)- Thus the above equation can be rewritten as 



v(x) X 



curl.xCurU: / g{y)G{y,x)ds{y) 



curUcurl^. / g{y)G{y,x)ds{y) 



+ 



x) X 



curl T cur Ir curl. 



+i\ 



curh^curU / g{y)G{y,x)ds{y) 



• f{x)ds{x) 
• f{x)ds{x) 
g{y)G{y,x)ds{y) 
} • h{x)ds{x) 



Making use of the vector identity: 

{curl^.curU.[5t(y)G(y,x)]} • h{x) = {curl j^curl j,[/i(x)G(y, x)]} ■ g{y), 
we obtain by a direct computation that for any g £ L^(Th), 



k- 



' [ E{y) ■ g{y)ds{y) = 0, 



where, for y G IK^\r, 



E[y) = -g^{cvii\ycm\y j G{y,x)f{x)xv{x)ds{x) — cm\ycm\y 



G{y,x)f{x) X i>{x)ds{x) 



+k curl. 



G{y,x)h{x) X i/{x)ds{x) + iXcuilyCurly / G{y,x)h{x)ds{x)} 



Thus we have 



curl curl E — k^E = 
vx E = {) 



It is clear that E{y) propagates upward above T satisfying the Rayleigh expansion radiation 
condition (jl.6p and propagates downward below F satisfying the Rayleigh expansion radiation 
condition (jl.6p with a replaced by —a. By the uniqueness of the radiating solution to the 
exterior problem of the Maxwell equations with the perfectly conducting condition and the 
analytic continuation of the solution of the Maxwell equations, it follows that E{y) = for 
2/3 > f{yi-,y2)- When y — t- F, the following jump relations hold on F: 



u X E^ — V X E 
iXE+ - iXE- 
V X curl E'^ — V X curl E~ 



0, on Fd, 
—iXh, on Tj, 
iXh, on F/, 



where the superscripts + and — indicate the limit obtained from $7 and M^\Q, respectively. It 
should be remarked that, since / G -ff~^/^(curl , F), the first integral over F in the definition of 
E{y) is well defined with //^^/^-density (see |21| ) and the corresponding jump conditions are 
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interpreted in the sense of the L^-hmit. Combining these jump relations and using the fact that 
V X = 1/ X curl E'^ = lead to 

curl curl S - A;^£' = 2/3 < 7(2/1,2/2), 

V X E'' =0 on T/), 

V X curli?" + iXE^ =0 on F/. 

A similar argument as in [15j can be applied to the above problem to show that E{y) = for 
2/3 < f{xi,X2)- Thus, we have 

/= [curlS]|r^ =0, h = -[iyx E]\r, =0, 

where [■]\rA stands for the jump across Ta of a function with A = D,I. The proof of Lemma 
I3.1l is thus completed. □ 



— 1/2 

The near field operator A'^ is defined by a bounded operator from B into (div, r^) 

which maps the boundary data (/ii,/i2) G -B to the tangential component 63 x E'^(x)\ti, 
near electric field. Here, E^ stands for the unique Rayleigh expansion radiating solution to the 
Maxwell equations with the following boundary conditions: 

1/ X E^ = hi on r_D, X curlE* - iA(ii;'*)T = /i2 on F/. 

By the well-posedness of the direct problem (see |15) ) it is known that is injective and bounded. 
Furthermore, A^ is a compact operator. To see this, we need the following periodic representation 
formula. 

Lemma 3.2. Assume that E satisfies the Rayleigh expansion radiation condition /il.6\) and the 
Maxwell equations in fl. Then for any a; G O we have 

E{x) = curU j i^{y) X E{y)G{x,y)ds{y) + ^curU.curU j v{y) x cm\ E{y)G{x,y)ds{y), 
where G{x,y) is the quasi-periodic Green function defined by Ii2.2\) . 

Proof. For arbitrarily fixed x G $7 and an arbitrary constant vector P G let F(x, y) = PG{x, y) 
with y G fi. Assume that x G 0^ for some 6 > 0. Denote by Bs(x) the small ball centered at 
X with radius 6 such that Bs{x) C Oft. It is clear that both E and F{x, •) satisfy the vector 
Helmhotz equation in Q\B^{x). Using Green's second vector theorem and the quasi-periodicity 
of E and F we have 



= [ E{y)-AF{x,y)-AE{y)-F{x,y)dy 
Jnt\Bs{x) 

= ("/ + /+ / ]{i^x E -curly F + i^-EdiVyF- uxF-curlE}ds{y) 

\ Jr JVu J\y-x\=5 



b ■J\y-x\=s^ 

h + h + h. (3.2) 
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By a direct computation, we have 



j {uxE- cml y[PG{x,y)] + v ■ Edivy [PG(x, y)] -vx [PG(x,y)] • curlE} ds{y) 
j {{i^xE)x VyG{x, y) + {y E)VyG{x, y) + {i^ x curl E)G{x, y)} ds{y) ■ P 

TG{E){y)ds{y) ■ P (3.3) 



with 

TG{E){y) = {ux E{y)) x VyG{x,y) + {u{y) ■ E{y))VyG{x,y) + {u x cml E{y))G{x,y) 

= curl,[i/(y) X E{y)G{x,y)] - VMv) ■ E{y)G{x,y)] + {u x cml E{y))G{x,y), 

where we have used the fact that \/yG{x,y) = —'VxG{x,y) to get the second equahty, 

TG{E){y)ds{y) ■ P 



= / {[{VyG{x,y)xP)xe3]-E{y) + iVyGix,y)-P){e3-E{y)) 

-PG{x, y) ■ (curl E{y) x eg) } ds{y) = 0, (3.4) 

where the last equality follows from the Rayleigh expansion condition (|1.6p for E{y), the definition 
of G{x, y) and the fact that divE = in 



h ■■= [ TG{E){y)ds{y) ■ P 

J \y—x\=& 

= it TG-^{E){y)ds{y) + f T^{E){y)ds{y)\ ■ P (3.5) 

\j\'ii—x\=5 J\'ii~x\=5 I 



Since G{x,y) — ^{x,y) is a C°°-function with respect to y in Bs{x) (see |23|). we have 
^\y_x\^^TG-^{E){y)ds{y) — )■ as 5 — )• 0. The application of the mean value theorem yields 
that f^y_^^^^Tq,{E){y)ds{y)) E{x) as J ^ (cf. [H]). Thus it follows from (I32|)-(l33|) that 



E{x) = cml^ j^v{y)x E{y)G{x,y)ds{y) + j^v{y)xcmlE{y)G{x,y)ds{y) 

-V,J^u{y)-Eiy)Gix,y)dsiy) 

which is analogous to the well-know non-periodic Stratton-Chu representation theorem ([1] 
Theorem 6.1]). Finally, the application of the Stokes theorem together with the vector identity 
curl curl = — A + V(V-) gives the desired result. □ 

It is seen from Lemma 13.21 and the well-posedness of the direct scattering problem that N is 
a composition of a bounded operator mapping the boundary data into the scattered field with a 
compact operator taking the scattered field to its tangential component of E^ on T;,. Thus is 
compact. We now prove, with the help of Lemma [3.21 that N has a dense range in (div, Ff,). 
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Lemma 3.3. The set {A^(v?,'0) I f e i^(rD),V' G LfC^l)} is dense in ^/^(div, Tf,). 
Proof. Let h G H^^^^ {cml ,rb) = i^^^^^div, Ffe)' satisfy 

<iv(v3,V'),/i>=o V(^Gy(rD), ^G^?(r/), (3.6) 

where < •,• > denotes the duahty between (curljFf,) and (div, Tb). Then it is 

sufficient to prove that h = 0. By the definition of N and the well-posedness of the direct 
scattering problem there exists a unique E € Hiodcml , fl) satisfying the Rayleigh expansion 
radiation condition (|1.6p such that N{(p,ip) = x E on Tb- From Lemma 13.21 it foUows that 



< N{{p,ilj),h > = I ez X E{x) ■ h{x)ds{x) 



{h{x) X 63) • jcurU j y{y) X E{y)G{x,y)ds{y) j ds{x) 
+^ j X 63) • jcurUcurU j u{y) x cuTlE{y)G{x,y)ds{y) \ ds{x) 



:= h + h- 
Interchanging the order of integration gives 



h = - J viy) X E{y) ■ jcurly j G{x,y)h{x) x esds{x) } ds{y), 
I2 = -j^ J I' (y) X curl E{y) ■ i^cuily curly J G{x,y)h{x) x e^ds^x) \ ds{y). 



Let 



F{y) := -^curlycurly f G{x,y)h{x) x esdsix), y G M^\Fb. 
Then, since curl curl = — A + V(V-), we have 

curlF(?/) = — curlj^ / h{x) x ej,G{x,y)ds{x). 

Thus 

< N{>f X ip),h >= / X E{y) ■ cml F{y) - u{y) x F{y) ■ cuvlE{y)ds{y). (3.7) 



Let E be the — a-quasi-periodic Rayleigh expansion radiating solution to the problem: 

curl curl ^ — A:^£' = in fi, 

X E = u X F on Fd, (3.8) 

X curl^ - iA^T = X curlF - iAFr on Fj. (3.9) 
From Green's second vector theorem and the Rayleigh expansion of E and E it follows that 

i/(y) X E{y) ■ cml E{y) - v{y) x E{y) ■ curl E{y)ds{y) 

v{y) X E{y) ■ curl E{y) - v{y) x E{y) ■ curl E{y)ds{y) = 0, 

b 
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which, in conjunction with the boundary conditions (|3.8p and (|3.9p . u x E = if on To and 
u X curli? — iXEt = ip on T/, imphes that 



u X F ■ curl E + {u X curl F - IXEt) ■ E = if- curl E+ ip ■ E. 



This together with (|3.7p yields 



< Niip,ij),h > 



/ ip ■ curl F — V X F ■ curl E + j ip ■ F — {v x curl F — iXEx) • E 
Jvn Jr, 



if ■ 



curl F — curl E 



+ I ij-{F-E). 



It is seen from the above identity that the conjugate operator of N is given by 



N*h = (^(curlF - curl^)T, {F - E)Tj G B. 



Combining (j3.6p and the boundary conditions ()3.8p and (|3.9p gives x F = x and x curl F = 
V X curl£^ on F. By Holmgren's uniqueness theorem, F = E\n 17;,. Since vxF = vxE onTi, and 
both F and E satisfy the — a-quasi-periodic Rayleigh expansion radiation condition for X3 > 5, 
it follows from the uniqueness result for the exterior Dirichlet problem that F = E for X3 > b. 
Now, in view of the fact that E is analytic in 0, we have by the jump relation of curlF(y) as 
2/ — 7- Ffe that 



h = [curli*""*" — curlF ] Ir^ = [curli?^ — curlE' ] Ir^, 
which completes the proof of the lemma. 



0, 



□ 



4 The Linear Sampling Method 

For g G L^(Ft) consider E^^(x;g) defined by (|3.ip as incident waves. Denote by E^{x;g) the 
scattered solution of the problem (ll.2p - (ll.5p corresponding to E^"'(x;g). To derive a periodic 
version of the linear sampling method consider the following near field equation: 

F{gz) := 63 X E^{x,gz)ds{x) = 63 X cuvl^cml^{PG{x,z)} on F^, (4.1) 

where z £ {z £ W^\0 < Z3 < b} and P € is a polarization vector. It is clear that 

NH{g) = -F{g). (4.2) 

Theorem 4.1. Assume that F is Lipschitz continuous with the dissection F = F/j U S U F/ and 
T/ / 0. 

(1) If z £ R.^\Q, then for any e > there exists a gl p £ (F^) such that 
WiJ^alp) -J^x curlcurl{PG(-,z)}||^2(r^) < e 

and 

llff2,pllL2(rt) ^ 00 as z^r-. 
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(2) If z G ^, then for any e > and 6 > there exists a g'f'p S L^(rfc) such that 
II Wp) - ^ X curlcurl{PG(-,z)}|L2(r^) < 6 + J 

and 

Mz%\\lI{v^)^ as (5^0. 

Proof. (1) Let z G M^\r2. In this case, 63 x curla;Curla;{PG(x, ^)}|rj, is in the range of since it 
is the tangential component of the electric field E\> ^ '■= curla;Curl2.{PG(x, z)} which is a solution 
of the exterior mixed boundary value problem with boundary data hi = v x Ep^ on F^) and 
/i2 = 1/ X curlii^p^ — i\{Ep .^)t on F/, that is, 

63 X curl curl {PG(x,z)}|r(, = iV(/ii,/i2)- (4.3) 

It can then be seen from the denseness of the range of H that, for every e>0 there is a gfp^ := 
g{-;e,P,z) £ Lti^b) such that 

\\Hig'p,z) + (^i,^2)||y(r^)xL?(r,) < e- (4-4) 
The boundedness of N implies that 

for some positive constant C. From this, (|4.2p and (|4.3p it follows that 

\mip) - V X curlcurl{PG(-,z)}||i2(r^) < Ce. 
Furthermore, if z — )• F~, then we have 

||(/ii,^2)||y(ro)xL2(r,) ^ 00 
due to the singularity of h\ and /i2 as z — t- F^. This, together with (|4.4p . gives rise to 

_^hm_ ||-H'5'p,J|y(rc)xL2(r,) = oo> 

which together with the boundedness of H implies that 

^lmi_ ||5p,JlL2{r6) = oo- 

(2) Let z G il. In this case, 63 x curl curl {PG(a;, z)}|ri, is not in the range of N since, 
otherwise, curl curl z)} will be a solution to the Maxwell equations in M^\ri which is 

impossible due to its singularity at z. However, using the Tikhonov regularization, we can con- 
struct a regularized solution to the near field equation (j4.ip since, by Lemmas 13.11 and 13. 3| J- 
is compact and has a dense range. Specifically, for an arbitrary 5 > there exist functions 
^\pz^^\.Pz) ^ ^(^d) X L^(F/) corresponding to some parameter a = a{5) chosen by a regu- 
larization strategy (e.g., the Morozov discrepancy principle) such that 

\\N{hip,,,hlp,,)-ux curlcurl{PG(-,z)}||^.(r,) < 5. (4.5) 
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Furthermore, using the regularization strategy and the Picard theorem (see |14| ) we get 

Yima{5) = 0, J™ ||(/ii,p,^, /i2,P,JI|y(rD)xL2{r,) = oo- (4-6) 

Then by Lemma [3. II and the boundedness of A^, for any e > it is possible to find a gp^^ G (F;,) 
such that 

\\N{Hg'fy - N{hip^,,hlp^MLUT,) < (4.7) 
Thus we have from (j4.5p and (j4.7p that 

ll(-^ffp') - " curl curl {PG(-,z)}|L2(r^) <e + 6. 

Finally, by (I4.6P and (I4.7P in conjunction with the boundedness of H and N, we have 

ll5^;pllL2(r,) ^ oo as 6-^0. 

The proof is thus completed. □ 

We now discuss some numerical strategies on the implementation of the above linear sampling 
method. 

As stated in Section [3l the incident waves E^^{x;g) defined by ()3.ip are not of physical 
relevance since they propagate away from the surface. Thus {x; g{y)), the scattered field 
corresponding to Eg"'{x), can not be generated directly. In what follows, we make use of the 
method of Arens and Kirsch [3j to generate E^{x;g). We first examine that 



G{x,y)-Giy,x) = ^| ^ -i-e^(°"-(^-2')-'^"(2'3-x3)) + ^ 1 gi(an-(x-y)+/3„(^^ 

:= A(^)(x;y)+A(^)(x;y) (4.8) 



for y £ Tf, and x £ ilb- Note that A('^)(x;y) and A^^\x;y) are upward and downward propa- 
gating modes respectively. Set 



E^^\x;g) := curl^curU. / g{y)A^^\x;y)ds{y), 

J Ft 

E(^)(x;5) := curl.curU / ^A(^)(x; y)(is(y), 
E'^ix;g) := curUcurl, ^ ^ {g(x; y) - A(^)(x, y)} ds(y). 

Clearly, E'^^{x;g) is propagating towards the scattering surface, so the corresponding unique 
scattered filed E^{x;g) can be computed directly. It is seen from (14. Sp and the boundary value 
of E'{x,g) that 

ux[E^ix;g)\r^+E(-^\x;g)\r^} = ,y x [s^ix; g)\r^ + E'^{x; g)\r^ - E'^ix; g)\r^} 

= -uxE'^^{x;g)\r,. 
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Similarly, we have 



V X curl 



{x-g) + (x; g)] - iX {x;g) + (x; g) 



V X curl 



E\x-g)^E^''{x-g)-E'^{x-g) - iX E'{x; g) + E'^ix; g) - E'^ix; g) 



= -{z^ X cml E'"" {x; g) -iXE'''{x;g)T} ■ 

It follows from the uniqueness of the direct scattering problem that E^{x;g) = E^{x;g) + 
E^^\x;g). Thus we can exactly generate E^{x,g) using the incident field E^'^{x;g). 

Remark 4.2. The reason why we use E^'^(x;g) defined by (|3.ip as incident waves to reconstruct 
the periodic structure is that, by Lemma [3. H such kinds of incident waves lead to a dense range 
in B of the operator H. It should be remarked that, if F/ = 0, that is, the total electric field 
E{x) satisfies the perfectly conducting boundary condition x = on F, then we are allowed 
to choose the following field as incident waves: 

:= curUcurU / g{y)G{x,y)ds{y), X3 < b 

with g G Lj(rb). This kind of incident fields leads to a dense range in H-^/^{dw,T) of the 
operator mapping g into the tangential component on T oi E^^{x]g). Since E^J^{x;g) propagates 
downward in fi^, the corresponding scattered field can be produced directly. Thus the above 
strategy of generating the scattered field corresponding to E^"'{x;g) can be avoided. 

Our reconstruction algorithm consists of the following three steps: 

Step 1. Select a mesh of sampling points in a computing region S;, = {x G K^jO < X3 < fe, < 
xi,X2 < 27r} which contains the grating surface. 

Step 2. Making use of the Tikhonov regularization and the Morozov discrepancy principle to 
compute an approximate solution g^ p to the near field equation (14. ID . 

Step 3. Consider 1151 pllia^r;,) ^ indicator function of the sampling points z and get the 
contour plot of \\gl pWi'^^Tt) ^ ^ function of z. 

Remark 4.3. The above algorithm has been implemented in [16j for the two-dimensional TE 
polarization case, where the Maxwell equations are replaced by the scalar Helmholtz equation and 
the boundary conditions on Fd and Tj are replaced with the Dirichlet and impedance conditions, 
respectively. The numerical reconstruction results presented in [16j have shown the efficiency of 
the algorithm. The implementation of the above algorithm for the three-dimensional case of the 
full Maxwell equations is still in progress. 

Remark 4.4. It follows form Theorem 14 . 1 1 that the indicator function can be used to characterize 
the different region below and above the grating surface. The numerical implementation of the 
linear sampling method can be found in [9J for inverse electromagnetic scattering problems by 
general bounded obstacles, which has been proven to be very successful and effective once the 
necessary direct scattering data are available. In order to get a better reconstruction result the 
mesh in Step 1 must be fine so that the characterization of the grating surface could be clear. 
But this also increases the computational cost since the near field equation (|4.ip must be solved 
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at each sampling point. To avoid this, a multilevel linear sampling method was proposed in |20| 
for the inverse acoustic obstacle scattering problems. The multilevel linear sampling method has 
been shown to be effective and to possess asymptotically optimal computational complexity and 
thus provides a fast numerical technique to implement the linear sampling method. 
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